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Abstract. Two different self-adjoint Pauli extensions describing a spin-1/2 two-dimen- 
sional quantum system with singular magnetic field are studied. An Aharonov-Casher type 
formula is proved for the maximal Pauli extension and it is also checked that this extension 
can be approximated by operators corresponding to more regular magnetic fields. 



1. Introduction 

Two-dimensional spin-1/2 quantum systems involving magnetic fields are described by the 
self-adjoint Pauli operator. One interesting question about such systems is the appearance of 
zero modes (eigenfunctions with eigenvalue zero). Aharonov and Casher proved in Pj that 
if the magnetic field is bounded and compactly supported, then zero modes can arise, and 
the number of zero modes is simply connected to the total flux of the magnetic field. Since 
then, Aharonov-Casher type formulas have been proved for more and more singular magnetic 
fields in different settings, see [BJEHElEil- Recently they were proved for measure- valued 
magnetic fields in [H] by Erdos and Vougalter. 

We are interested in the Pauli operator when the magnetic field consists of a regular part 
with compact support and a singular part with a finite number of Aharonov-Bohm (AB) 
solenoids 0. The Pauli operator for such singular magnetic fields, defined initially on smooth 
functions with support not touching the singularities, is not essentially self-adjoint. Thus there 
are several ways of defining the self-adjoint Pauli extension, depending on what boundary 
conditions one sets at the AB solenoids, see P G3 ^2 EI- Different extensions describe 
different physics, and there is a discussion going on about which extensions describe the real 
physical situation. 

There are two possible approaches to making the choice of the extension: trying to describe 
boundary conditions at the singularities by means of modelling actual interaction of the particle 
with an AB flux, or considering approximations of singular fields by regular ones, see [51 118j . 
We are going to study the maximal extension introduced in called the Maximal Pauli 
operator, and compare it with the extension defined in 8j, that we will call the EV Pauli 
operator. These two extensions were recently studied in |16| in the presence of infinite number 
of AB solenoids, and it was proved that a magnetic field with infinite flux gives an infinite- 
dimensional space of zero modes for both extensions. 

When studying the Pauli operator in the presence of AB solenoids one must always keep 
in mind the possibility to reduce the intensities of solenoids by arbitrary integers by means of 
singular gauge transformations. In Section [21 we define both extensions via quadratic forms. 
The Maximal Pauli operator can be defined directly for arbitrary strength of the AB fluxes, 
while the EV Pauli operator has to be defined via gauge transformations if the AB intensities 
do not belong to the interval (—1, 1). The EV Pauli operator is not gauge invariant. However, 
following we always make a reduction of the AB intensities to the interval [—1/2,1/2). 
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Hence the EV Pauli operator is not uniquely denned for AB intensities belonging to (— 1, 1) \ 
[—1/2, 1/2), see Sectional Moreover, the asymmetry of the interval [—1/2, 1/2), leads to the 
absence of the invariance of spectral properties of the EV Pauli operator under the changing 
sign of the magnetic field, the latter property being natural to expect. 

For the Dirac operators with strongly singular magnetic field the question on the number of 
zero modes was considered in The definition of the self-adjoint operator considered there 
is close to the one in Erdos-Vougalter, however it is not gauge invariant, therefore the Aharonov 
Casher-type formula obtained in ^3] depends on intensity of each AB solenoid separtely. 

In Section [21 we establish that the Maximal Pauli operator is gauge invariant and that 
changing the sign of the magnetic field leads to anti-unitarily equivalence. Our main result is 
the Aharonov-Casher type formula for the Maximal Pauli operator. An interesting fact is that 
this operator can have both spin-up and spin-down zero modes, in contrary to the EV Pauli 
operator and the Pauli operator for less singular magnetic fields, which have either spin-up or 
spin-down zero modes, but not both. In (lOj a setting with an infinite lattice of AB solenoids 
with equal AB flux at each solenoid is studied, having both spin-up and spin-down zero modes, 
both with infinite multiplicity. 

In Section 01 we discuss the approximation by more regular fields in the sense of Borg and 
Pule, see j^J- It turns out that the Maximal Pauli operator can and the EV Pauli operator 
can not be approximated in this way. However, different ways of approximating the magnetic 
field may lead to different results, see 0] EH . 

2. Definition of the Pauli operators 
The Pauli operator is formally defined as 

P = (a ■ (-<V + A)) 2 = (-*V + A) 2 + a 3 B 
on L2(K 2 ) <g> C 2 . Here a = (a±, 02), where ci, a% and 03 are the Pauli matrices 

VI = (j J) , °2= ( • 0*) ' ^ ^ = ((J -l) ' 

A is the real magnetic vector potential and B = curl(A) is the magnetic field. This definition 
does not work if the magnetic field B is too singular, see the discussion in |SJE1- If A £ 
i2,ioc(R 2 ), using the notations = — idt + Ak, for k = 1,2, Q± = Hi ± HI2 and A for the 
Lebesgue measure, the Pauli operator can be defined via the quadratic form 

(1) ?#]HIQ + V + I! 2 + I!Q-?M 2 = fw-(-iV + A)Tp\ 2 d\(x), 

the domain being the closure in the sense of the metrics p[ip] of the core consisting of smooth 
compactly supported functions. With this notation, we can write the Pauli operator P as 



(2) P 









0+ 


I) 





Q\Q+ 
Q*_Q- 



However, defining the Pauli operator via the quadratic form p[ip] in Q requires that the 
vector potential A belongs to L2.io C (K 2 ), otherwise p[%li\ can be infinite for nice functions ip, 
see ■ If t ne magnetic field consists of only one AB solenoid located at the origin with 
intensity (flux divided by 2tt) a, then the magnetic vector potential A is given by A(xi, X2) = 
x i" x i (~^2 ; Xi) which is not in i2,ioc(K 2 ). Here, and elsewhere we identify a point (£1,2:2) in 
the two-dimensional space IR 2 with z — x\ + 1x2 in the complex plan C. 
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Following 0, we will define the Pauli operator via another quadratic form, that agrees with 
p[ip] for less singular magnetic fields. We start by describing the magnetic field. 

Even though the Pauli operator can be defined for more general magnetic fields, in order 
to demonstrate the main features of the study, without extra technicalities, we restrict ourself 
to a magnetic field consisting of a sum of two parts, the first being a smooth function with 
compact support, the second consisting of finitely many AB solenoids. Let A = be a 

set of distinct points in C and let otj £E \ Z. The magnetic field we will study in this paper 
has the form 



(3) B(z) = B (z) +J22™ j 5 Zjt 

j=i 

where Bq 6 Cq(R 2 ). In |S] the magnetic field is given by a signed real regular Borel measure 
/i on R 2 with locally finite total variation. It is clear that (x = Bo(z)dX(z) + Sj=i 27rOj<5 Zj is 
such a measure. 

The function h given by 



satisfies A/i = B since B e Cq(R 2 ) and A log \z — Zj\ = 2ir8 Zj in the sense of distributions. 
The function 

n 

h(z) = h (z) + ^2 a j l°g \ z ~ z ]\ 

i=i 

satisfies Ah = B. It is easily seen that h (z) ~ $ log \z\ as \z\ —> oo, and thus the asymptotics 
of e h W is 



e ±h(z) 




\z\ —> oo 

±a 3 



z — z,,\ 3 , z — > z 
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where $ = ± J B (z)d\(z) and $ = i f B{z)d\(z) = $ + J]™ =1 a r 

We are now ready to define the two self-adjoint Pauli operators. The decisive difference 
between them is the sense in which we are taking derivatives. This leads to different domains, 
and, as we will see in later sections, to different properties of the operators. Let us introduce 
notations for taking derivatives on the different spaces of distributions. Remember that A = 
{zj}" =1 is a finite set of distinct points in C. We let the derivatives in P'(R 2 ) be denoted by 
d and the derivatives in V(M. 2 \ A) be denoted by d with a tilde over it, that is d. Thus, for 
example, by d z we mean in the space Z?'(K 2 ) and by d z we mean in the space X>'(R 2 \ A). 

2.1. The EV Pauli operator. We follow [B] and define the sesquilinear forms ir + and 7r_ by 

- 4 J 8- z (e-V+R (e- fe e+) e 2h dX(z), 
D(tt^) = G L 2 (R 2 ) : 7r^(^+,V+)<oo}, 

and 



B(i h ) = {f £l 2 (M 2 ) : 7r?L(V_,V_) < oo} 
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Now set 



D(tt' 1 ) =2)(^)©D(^) = {V> = fjM £l 2 (t 2 )®C 2 : tt^^Oc}- 

Let us make more accurate the description of the domains of the forms ir± and ir h . For 
example, what is required of a function ip + to be in D(tt+)? It should belong to ^(K 2 ), and 
the expression 

^ (</>+>+) =4 [\d- z (e~ h ^ + )\ 2 e 2h d\(z) 



should have a meaning and be finite. This means that the distribution d% (e h ip+) actu- 
ally must be a function and its modulus multiplied with e h must belong to L2(K 2 ), that is 
\dg {e~ h ip+) \e h G L2(M 2 ). This forces all the intensities ctj to be in the interval (—1, 1), see 
Next we define the norm by 

HivC*. = \H+\\\lx + 111^-111^, 

where 

\U + \\\l h+ = H + \\ 2 + \\d- z (e- h ^ + )e h \\ 2 

and 

= ||^|| 2 + ||a z (eV-)e-'i| 2 . 

This form TT h is symmetric, nonnegative and closed with respect to || • ||, again see [H], and 
hence it defines a unique self-adjoint operator Vh via 

(4) I>CP ft ) = U> G T>{n h ) : tt^, •) G (L 2 (R 2 ) ® C 2 )} 
and 

(5) (^,0=^,0. ^ e D(7\U g 2)(tt' 1 ). 

We call this operator Vh the non-reduced EV Pauli operator. 

If some intensities ctj belongs to R \ [—1/2, 1/2), we let a* be the unique real number in 
[—1/2, 1/2) such that ctj and a* differ only by an integer, that is a* — ctj = m,j G Z. We define 
the reduced EV Pauli operator (or just the EV Pauli operator), Ph, to be 

(6) P h = exp(i(j))Vh exp(-icj)) 

where <p{z) — X)j l =i m j ar g( z — z j)- Hence, if there are some otj outside the interval (—1, 1) 
only the reduced EV Pauli operator is well-defined. If all the intensities ctj belong to the 
interval [—1/2, 1/2) then we do not have to perform the reduction and hence there is only 
one definition. However, if there are intensities ctj inside the interval (—1, 1) but outside the 
interval [—1/2, 1/2) then we have two different definitions of the EV Pauli operator, the direct 
one and the one obtained by reduction. In the next section we will show that these two 
operators are not the same. 
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2.2. The Maximal Pauli operator. Now, again, let aj G R \ Z. We define the forms 



and 



4 J d- z (e-' l e+) e 2h d\(z), 

p h _W-,Z-) = 4 J B z (ehp7)d z (e h ^) e- 2h d\(z), 
V(p h _) = G L 2 (R 2 ) : p^(V-,V-) < oo}. 



Now set 



D(p /l ) 



D(p^©D(p1) = {^= ( J+j GL 2 (K 2 )«C 2 : p"^) < oc}. 



Again, let us make clear about the domains of the forms. For a function ip + to be in T>(p^ ) it 
is required that ip + G L 2 (R 2 ) and that d z {e~ h ip + ) is a function. After taking the modulus of 
this derivative and multiplying by e h we should get into L 2 (K 2 \ A), that is \d z {e~ h ip + )\e h G 
L 2 (R 2 \ A). Note that the form p h does not feel the AB fluxes at A since the derivatives are 
taken in the space V'(M. 2 \ A), and integration does not feel A cither since A has Lebesgue 
measure zero. This enable the AB solenoids to have intensities that lies outside (—1, 1). 
Also, define the norm 



111^ 



where 



and 



IIIVm 



d- z (e"V+) e h 



H-\\ 2 + ~d z {e h i,_). 



Proposition 2.1. The form p h defined above is symmetric, nonnegative and closed with respect 
to || • ||. 

Proof. It is clear that p h is symmetric and nonnegative. Let ip n = (ip n ,+ ,ipn,-) be a Cauchy 
sequence in the norm |||-||| p h. This implies that ip n ,± — » V'i m L 2 (d\(z)), d z (e~ h ip n , + ) — > u + in 
L 2 (e 2h d\{z)) and d z (e h ip n -) -> u_ in L 2 {e~ 2h d\(z)) . We have to show that d 2 (e _h ^+) = «+ 
and d z (e h ip-) = u— For any test-function G C^R 2 \ A), 



1 0(u+-a s (e-V+))dA(«) 



< j 4>(u+- d- z (e- 
+ | 4(<A)e"' 1 (V+ - Vn,+) 

< ||^1 • u+ - a 2 (e-Vn,+) 



L 2 ( e 2h) 



+ 



9z^)e-" •||V + -^+l|. 
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The last expression tends to zero as n — » oo, since the first terms in each sum is bounded 
(thanks to 0) and the other one tends to zero. The proof is the same for the spin down 
component. This shows that p h is closed. □ 

Hence p h defines a unique self-adjoint operator tyh via 

(7) D(«p h ) = {^D(p' 1 ) : p' l (W) G (L 2 (K 2 )®C 2 )} 
and 

(8) Whil>,Z)=P h y>,t), ^^Ue^)- 

We call this operator tyh the Maximal Pauli operator. 

3. Properties of the Pauli operators 

In this section we will compare some properties of the two Pauli operators Ph and tyh 
defined in the previous section. We start by showing that tyh is gauge invariant while Ph is 
not. 

3.1. Gauge transformations. Let B(z) = P>q(z) + Y^ = i ^(Xj5 Zj be the same magnetic field 
as before and let B{z) be another magnetic field that differs from B(z) only by some multiples 
of the delta functions, that is B(z) — B{z) = Xw=i 27rTO,-<5 Z;j , where rrij are integers, not all 
zero. Then the corresponding scalar potentials h(z) and h(z) differ only by the corresponding 
logarithms h(z) — h{z) — 5Zj=i m j 1°S \ z ~ z j\- Now with 4>{z) — Y^j=i m j ar g( 2 ~ z j) we § e t 
h(z) + i(j)(z) = h(z) + Y)j—i nij \og(z — zj). This function is multivalued, however, since rrij 
are integers, we have 



n 

(9) d E (h(z) + i<f>(zf) = dzh(z) + mjd E log(z - zj), 

(10) dg (h(z) + i<t>(z)) = d E h(z), and 



(11) e h+i * = e h ]J(z-z j ) m K 

To see that Ph is not gauge invariant it is enough to look at an example. Let n = 1, Z\ = 0, 
a,\ = —1/2 and m\ — 1, so the two magnetic fields are B(z) — B (z) — it5o and B(z) = B (z) + 
ttSo. The scalar potentials are given by h(z) = h (z) — ~log\z\ and h(z) = h (z) + ^\og\z\ 
respectively, where h (z) is a smooth function with asymptotics <& log|z| as \z\ — » oo. We 
should show that T)(n h ) is not given by e- t4 "D(TT h ) 7 where <f>{z) = arg(z). Then it follows that 
TT h and ir h do not define unitarily equivalent operators. 

Let e D(jt}). This means, in particular, that dz(i/j+e h ) belongs to Li i i oc (R 2 ). Now 
let ip+ = e~^ip+. Then, according to 1)11(1 we get 

dS+e-' h ) = dnty+e-*-*) = d z {^—^j = W+e~^ + t/> + e- h v6 

which is not in Li y \ oc (M. 2 ) since it is a distribution involving So (for non-smooth it is not 
even well-defined). Thus we have D(tt^) ^ e'^V^) and hence V(w ft ) ^ e~ ll >"D{-K h ) so ir h 
and ir h are not defining unitarily equivalent operators. 

Let us now study what happens with p h when we do gauge transforms. Let ip = (ip+, "0-)' S 
D{p h ). We should check that e - ^V belongs to T)(p h ), where </>(z) = £™ =1 rrij arg(z — Zj) is the 



ZERO MODES OF THE PAULI OPERATOR 



7 



harmonic conjugate to h(z) — h(z). We do this for pif. It is similar for p'f. Since ip + 6 £>(p+) 

we know that dz(ip+e-' 1 ) E L Uoc (R 2 \ A). Let us check that 9j(^ + e"' 1 ) € L Uoc (]R 2 \ A). 
Again, by (|11J) we have 



3 = 1 



which clearly belongs to Li ! i oc (R 2 \ A). 

Next we should check that <9z(-0+e _ft, )|e' t belongs to L2(M. 2 \ A) under the assumption that 
\dz(ip + e~ h )\e h belongs to L 2 (IR 2 \ A). A calculation using l(TU|l and l(TT|l gives 



8- z e->+ e = \dM e 



(12) 



0,(-A(*)W>+ + *W+(*) 



3=1 



eh n i* ~~ 



|d 2 h i>+^ | e'\ 



Hence V>+ e TJ(p+) implies = e -i<! V+ e 2)(p+). In a similar way it follows that ip- G X>(p^_) 
implies that ip_ = e"^V- G ©(P-)- Thus e-^D(p h ) C D(p' 1 )- In the same way we can 
show that e^D^) C D(p' 1 ), and thus we can conclude that e -^2)(p' 1 ) = D(p' h ). From the 
calculation in 112|) and a similar calculation for -0_ it also follows that 



9, 



5 



„2/i 



,-2h 



dA(z) 



V-4 



d z (eV-) e~ 2/l dA(z) 



p*(lM)- 



Hence we can conclude that if tp e D(^) and £ 6 D(p' 1 ) then e"^ e and e ~^£ G 

CD(p' 1 ). If we denote by TJ$ the unitary operator of multiplication by e^, then we get 

and hence tyh and are unitarily equivalent. We have proved the following proposition. 

Proposition 3.1. Let B and B be two singular magnetic fields as in (j3J, with difference 
B — B = Y^j=i 2TrrrijS Z:j , where rrij are integers, not all equal to zero. Then their corresponding 
Maximal Pauli operators defined by Q and (JSJ are unitarily equivalent. 

3.2. Zero modes. When studying spectral properties of the operator tyh it is sufficient to 
consider AB intensities ctj that belong to the interval (0, 1), since the operator is gauge invari- 
ant. See the discussion after the proof of Theorem 13.31 for more details about what happens 
when we do gauge transformations. 
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Lemma 3.2. Let Cj G C and Zj G C, j = 1, . . . , n, where Zj ^ zi if j ^ % and not all Cj are 
equal to zero. Then 



(13) E 



z 



where I is the smallest nonnegative integer such that Y^j=i c j z j 7^ 0- 
Proof. If z| is large in comparison with all \zj\ we have 




0(\z\- 1 - 2 ) 

and thus Yl*j=i T=z~ ~ I 21 ! - ' -1 as M ~^ 00 • ^ 

Remark. We note that Z in Lemma [3.21 may never be greater than n — 1. Indeed, if Z > n 
then we would have the linear system of equations {X}j=i ^'^y = 0}fc=o- But the determinant 
of this system is fli^C^ — z j) 0, and this would force all Cj to be zero. 



Note also that for Z < n we have a system of Z equations {X^=i Cj 2 ? = 0}fc=o w ith 



unknowns Cj, and that the Z x n matrix {^} has rank Z. ■ 

Theorem 3.3. Let B(z) be the magnetic field © with all aj G (0, 1), and let tyh be the Pauli 
operator defined by J7J) and (jSJ in Section^ corresponding to B{z). Then 

(14) dim kcr <p h = {n -$} + {$}, 

where $ = J £?(z)dA(z), and {x} denotes the largest integer strictly less than x if x > 1 
and i/x < 1. Using the notations Q± introduced in Section\^ we also have 

(15) dimker<2+ = {n ~ $} and dimker Q_ = {<&} . 

Proof. We follow the reasoning originating in with necessary modifications. First we note 
that tp-Y belongs to ker^ if and only if tp + belongs to kerQ + and ip- belongs to kerQ_, 
which is equivalent to 

(16) Bs (e _ V+) = and d z (e h V-) = 0. 

This means exactly that f+(z) = e~ h t/j + (z) is holomorphic and f-(z) — e h ip_(z) is antiholo- 
morphic in z G M 2 \ A. It is the change in the domain where the functions are holomorphic 
that influences the result. 

Let us start with the spin- up component rp + . The function /+ is allowed to have poles of 
order at most one at Zj, j = 1, ...,n, and no others, since e h ~ \z — Zj\ aj as z — ► Zj and 
= f+e should belong to L2(R 2 ). Hence there exist constants Cj such that the function 
f+(z) — X)j=i T^7~ IS en ti re - From the asymptotics e h ~ |z|*, \z\ — > oo, it follows that 

— Y^j=i z-z- ma y om y be a polynomial of degree at most N = — $ — 2. Hence 



f+(z) = h a + aiz + ...(/ \ : ' 

Z — Zj 
3=1 J 
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where we let the polynomial part disappear if N < 0. Now, the asymptotics for tp+ is 

^ + (z)~izr i - i+ * + izr+*, i*i->oo, 

where I is the smallest nonnegative integer such that Y^j=i c j z j 0- To have ip+ in L2(K 2 ) 
we take I to be the smallest nonnegative integer strictly greater than $. Remember also 
from the remark after Lemma [3.21 that I < n — 1. We get three cases. If $ < — 1, then all 
complex numbers Cj can be chosen freely, and a polynomial of degree {— $} — 1 may be added 
which results {n — $} degrees of freedom. If— l<<f><n— lwe have no contribution from 
the polynomial, and we now have to chose the coefficients Cj such that Y^j=i c j z j = f° r 
k = 0, 1, . . . , I — 1. The dimension of the null-space of the matrix {zj} is n — I = {n — <f>}. If 
<& > n — 1 then we must have all coefficients Cj equal to zero and we get no contribution from 
the polynomial. Hence, in all three cases we have {n — $} spin-up zero modes. 

Let us now focus on the spin-down component The function /_ may not have any 
singularities, since the asymptotics of e~ h is \z — Zj\~ aj as z — > Zj. Hence /_ must be entire. 
Moreover, /_ may grow no faster than a polynomial of degree $ — 1 for -0_ to be in ^(K 2 ). 
Thus /_ has to be a polynomial of degree at most {$} — 1, which gives us {$} spin-down zero 
modes. □ 

The number of zero modes for tyh and Ph are not the same. The Aharonov-Casher theorem 
for the EV Pauli operator (Theorem 3.1 in jSj) states for the field under consideration: 

Theorem 3.4. Let B(z) be as in @ and let B(z) be the unique magnetic field where all AB 
intensities ctj are reduced to the interval [—1/2, 1/2), that is B(z) = B{z) + X)j=i ^ 7rm j^z j , 
where ctj + rrij 6 [—1/2, 1/2). Let <I> = j- J B{z)d\(z) . Then the dimension of the kernel of 
the EV Pauli operator Ph is given by {|$|}. All zero modes belong only to the spin-up or only 
to the spin-down component (depending on the sign of 

Below we explain by some concrete examples how the spectral properties of the two Pauli 
operators tyh and Ph differ. 

Example 3.5. Since Ph is not gauge invariant we must not expect that the number of zero 
modes of Ph is invariant under gauge transforms. To see that this property in fact can fail, let 
us look at the Pauli operators Ph ± and Ph 2 induced by the magnetic fields 

B\{z) = Bq(z) + ttSq, and 
B 2 (z) = B Q {z) - ttS 

respectively, where Bq has compact support and <I>o = ^- J Bo(z)dX(z) = |. Then Bi is 
reduced (that is, its AB intensity belong to [—1/2, 1/2)) but B\ has to be reduced. Due to 
Theorem 13.41 the EV Pauli operators Pf ll and Ph 2 corresponding to B\ and B 2 have no zero 
modes. However, a direct computation for the non-reduced EV Pauli operator Vh x correspond- 
ing to B\ shows that it actually has one zero mode. The situation is getting more interesting 
when we look at the operator that should correspond to B 3 = B {z) + 3-7r<5o. The AB intensity 
for B 3 is too strong so we have to make a reduction. In [5] the reduction is made to the interval 
[—1/2, 1/2), and we have followed this conventions, but physically there is nothing that says 
that this is the natural choice. Reducing the AB intensity of B% to —1/2 gives an operator 
with no zero modes and reducing it to 1/2 gives an operator with one zero mode. 

The Maximal Pauli operators ?phi> *P/i 2 an d Vh 3 for these three magnetic fields all have 
one zero mode. This is easily seen by applying Theorem 13.31 to tyht arid then using the fact 
that the operators are unitarily equivalent. 
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However, more understanding is achieved when looking more closely at how the eigenfunc- 
tions for these three Maximal Pauli operators look like. Let hk be the scalar potential for Bk, 
k = 1,2, 3. Then, as we have seen before hi(z) = ho(z) + \ log \z\, h 2 (z) = h (z) — \ log \z\ and 
h%(z) = ho(z) + I log \z\ where ho(z) corresponds to B (z). Following the reasoning from the 
proof of Theorem 13. 31 we see that the solution space to tyhxip — is spanned by ip = (0, e~' 11 )*. 

Next, we see what the solutions to *}3/ l2 tp = look like. Now we have <J> 2 = 57 / B 2 (z)dX(z) = 
1/4 > 0. Let us begin with the spin- up component ip+. This time, the holomorphic /+ = 
e~ h2 ip + may not have any poles since then ip + would not belong to L 2 (M 2 ), and f+(z) = 
e~ h2 ip + (z) — > as \z\ — > 00, so we must have /+ = 0, and thus ip + = 0. For ip-(z) to 
be in L 2 (M 2 ) it is possible for /_ to have a pole of order 1 at the origin. Hence there ex- 
ist a constant c such that f—(z) — c/z is antiholomorphic in the whole plane. The function 
f-(z) — > as \z\ — > 00 since the total intensity <f> 2 > 0. This implies, by Liouville's theorem, 
that f-(z) = c/z, so the solution space to tyh 2 *P = is spanned by ip{z) = (0, e~ h2 /z). 

Finally, let us determine the solutions to ySh 3 4> = 0- Now $3 = ^ J B%{z)d\(z) = 9/4. 
Consider the spin- up part ip+. For ip + to be in ^(K 2 ) our function / + may have a pole of 
order no more than two at the origin. As before, there exist constants c\ and c 2 such that 
f+{z) —ci/z — C2/z 2 is entire and its limit is zero as \z\ — > 00, and thus f+(z) = c±/z + 02/ z 2 . 
Again, both ci and c 2 must vanish for ip+ to be in L 2 (R 2 ) (otherwise we would not stay in L 2 
at infinity). Thus ip+ = 0. On the other hand, the function /_ may not have any poles (these 
poles would push ?/>_ out of L 2 (IR 2 )), so it is antiholomorphic in the whole plane. It also may 
grow no faster than |z| 5 / 4 as \z\ — ► 00, and thus /_ has to be a first order polynomial in z, 
that is f-(z) = c + C\Z. Moreover for ip_ to be in L 2 (M 2 ) it must have a zero of order 1 at 
the origin, and thus f-(z) = C\z. We conclude that the solutions to $Pft 3 V> = are spanned by 
(0,ze^' l 3)*. " ■ 



A natural property one should expect of a reasonably defined Pauli operator is that its spec- 
tral properties are invariant under the reversing the direction of the magnetic field: B — B. 
The corresponding operators are formally anti-unitary equivalent under the transformation 
ip 1 — ► ip and interchanging of ip+ and ip— . 



Example 3.6. The number of zero modes for P/, is not invariant under B(z) —B(z), 
which we should not expect since the interval [—1/2, 1/2) is not symmetric. We check this by 
showing that the number of zero modes are not the same. To see this, let B{z) = Bq(z) + tt5o, 
where Bo has compact support and ^0 = ^1 Bo(z)dA(z) = |. Then B has to be reduced 
since the AB intensity at zero is 1/2 ^ [—1/2, 1/2). After reduction we get the magnetic field 
B(z) = Bo(z)-nS , and we can apply TheoremEU Now let $ = ^ / B*d\(z) = \. Thus the 
number of zero modes for Ph is 0. Now look at the Pauli operator defined by the magnetic 
field B-(z) — —B(z) — ~Bq(z) — irdo. This magnetic field is reduced and thus we can apply 
Theorem 13.41 directly. The total intensity is $_ — ^ J —B{z)d\{z) = — |, so the number of 
zero modes for P—h is 1. If B has several AB fluxes then the difference in the number of zero 
modes of Ph and P-h can be made arbitrarily large. 

Now, let us check that the number of zero modes for tyh is invariant under B{z) 1— > —B(z). 
Since it is clear that the number of zero modes is invariant under z 1— > z we look instead 
at how the Pauli operators change when we do B(z) 1— > B(z) = —B(z). If we set ( = z 
we get B(C) = —B{z) and the scalar potentials satisfy h(Q = —h(z). Now assume that 
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%j> = (ijj + (z),iP^{z)) t G T)(p h ). Then 

( (i>+i z )\ (i>+( z ) 



= 4^ |9 z -(V+( 2 )e-' l(z) )| 2 c 2 ' l(2) + |9 z (V-(2)e' l(z) | 2 e" 2 ' l<z) dA(z) 
= 4 1 |a e (^ + (C)e^«| B e- a ^« + |a f (^_(C)e- Rra | a e 2t(C) dA(C) 

Hence we see that (■^ + ,^_) < belongs to !D($Pk(>)) if and only if (^_,^> + )* belongs to f^P/^) 
and then = ^ h(z) V where V : L 2 (R 2 ) ® C 2 -> L 2 (M. 2 ) <8> C 2 is the isometric operator 

given by V((^>+, tp-Y) = (V'-jVM-)'- Hence it is clear that tyu^ and ^Ph(z) have the same 
number of zero modes. ■ 

Example 3.7. In the previous example we saw that the number of zero modes for the Maximal 
Pauli operators corresponding to B and —B are the same. This can easily be seen directly 
from the Aharonov-Casher formula in Theorem 13.31 To be able to apply the theorem to 
— B = —Bo — Y^j=x ^ctjSj we have to do gauge transformations, adding 1 to all the AB 
intensities, resulting in B = —Bo + Y^j=i 2tt(1 — otj)Sj. Now according to Theorem 13.31 the 
number of zero modes of ty-h is equal to 

dimker*p_/j = {$} + {n - 4} = {n - $} + {$} = dimker^, 

where we have used that $ = j- J BdX(z) = n — $. ■ 



4. Approximation by regular fields 

We have mentioned that the different Pauli extensions depend on which boundary conditions 
are induced at the AB fluxes. Let us now make this more precise. Since the self-adjoint 
extension only depends on the boundary condition at the AB solenoids it is enough to study 
the case of one such solenoid and no smooth field. For simplicity, let the solenoid be located 
at the origin, with intensity a £ (0, 1), that is, let the magnetic field be given by B — 2naSo- 
We consider self-adjoint extensions of the Pauli operator P that can be written in the form 



P = 



P+ \ (Q* + Q+ 
P-J I QLQ-, 



with some explicitly chosen closed operators Q±. It is exactly such extensions P that can be 
defined by the quadratic form JIJ. A function tp + belongs to D(P + ) if and only if 4>+ belongs 
to D(Q+) and Q+ip+ belongs to D(Q* + ), and similarly for P_. 

With each self-adjoint extension P± = Q*±Q± one can associate (see 03 ^2 E|) func- 
tional ct a , c±, c^_ x and cf_ a , by 



12 



M. PERSSON 



c± a (V±) = lira ;r°l / ^±40, 

c±(^±) = Iimr- Q f ^ / - r- a ct(i>±) 

1 / >27r 

c^_i(^±) = limr 1 -"— / 7/>±e lf, d0, and 

r— »0 Z7T Jq 



such that ^)-t- G D(P±) if and only if 

(17) V± ~ c± a r- Q + c±r Q + c^r^e^ + c^r 1 "^"^ + 0{f) 

as r — > 0, where 7 = min(l + a, 2 — a) and z = re 10 . 

Any two nontrivial independent linear relations between these functionals determine a self- 
adjoint extension. In order that the operator be rotation-invariant, none of these relations may 
involve both a and 1 — a terms simultaneously. Accordingly, the parameters v = / c^ a 
and v± = cf-a/Ca-ii with possible values in (—00, 00], are introduced in [5], and it is proved 
that the operators P± can be approximated by operators with regularized magnetic fields in 
the norm resolvent sense if and only if i/^ = 00 and G (—00, 00) or if Vq G (— 00, 00) and 
= 00. We are now going to check what parameters the Maximal and EV Pauli operators 
corresponds to. 

Generally, for the function ip+ to be in D(P + ), it must belong to D(Q + ) and Q+ip+ must 
belong to D(Q!J_). We will find out what is required for a function g to be in D(Q1). Take 
any (f> + G D(Q+), then the integration by parts on the domain e < \z\ gives 



(.9,Q + 0+> = lim / g(z)-^(e-^ + (z))e h d\(z) 

£ ^°J\z\>e <JZ 

= - lim / |-( 3 (z)e' i )e-' l ^)dA(z) 

J\z\>e oz 

F /-27T 

-lim-/ g(se l0 )(j) + (ee w )e- l0 de 



= (-0-5,0+) - lim § g{ee l6 )<j> + {ee^)e- l0 d6 

e^O 2 ,/ 

Hence, for g to belong to D(Q1) it is necessary and sufficient that 

lim £ / g(ee w )(j) + {ee i9 )e- l9 d6 = 
for all </> + G D(p + ), and thus for Q+ip+ to belong to T)(Q* + ) it is necessary and sufficient that 

lime/ ( ^-( e - h ^ + )e h \ <i>+{se e )e- ie d6 = 

for all </>+ G D(p+). We know that ■0+ has asymptotics tp+ ~ cl Q 7-~ Q + c+7 ,Q + c+_ 1 r Q ~ 1 e~ je + 
c~l_ a r 1 ~ a e~ l6 + 0(r 7 ) and thatj| = -I- ~ in polar coordinates. A calculation gives 
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e— (e- h i>+)e h e 
oz 

hence we must have 



iO 



-2ac + _ a e~ a + 2(1 - a)c+_ Q e 1 - Q e 



(18) lim / (-2ac± a e- a + 2(1 - ay+^e^e- 19 ) <j>+(ee i9 )d9 = 
for all <p + G D(p + ). A similar calculation for the spin-down component yields 

(19) lim / (2ac-e a + 2(a~l)c-_ l e a ' 1 e %e )(t>-{ee ie )de = {). 

We will now calculate what parameters and uf 1 the Maximal and EV Pauli extensions 
correspond to. To do so, it is enough to study the asymptoics of the functions in the form 
core. 

Let us first consider the Maximal Pauli extension. Functions on the form (cf>Qc/z)e h con- 
stitute a form core for p , , where 4>o is smooth. Hence there are elements in D(p + ) that 
asymptotically behave as r a and also elements with asymptotics r a ~~ 1 e~~ l6 . According to (|18fl 
this means that c_ a and c^"_ Q must be zero. Similarly the elements that behave like r~ a and el- 
ements that behave like r 1 ~ a e ie constitute a form core for p'i , which by l|19fl forces c~ and c~_ ± 
to be zero. The parameters u and are given by z/q~ = c^/c~t a = oo, = c\_ a j c\_ Y = 0, 
isq = c~/cZ a = and v± = c^_ a /c~_ 1 = oo. Hence the Maximal Pauli operator can be 
approximated in the sense of [JJ. 

Let us now consider the EV Pauli extension, and study the case when a <E (0, 1/2). The 
case a < follows in a a similar way. A form core for 7r^ is given by e h 4>o where <j>o is smooth, 
see [H] ■ These functions have asymptotic behavior r a . From IjlSjl follows that ct a must vanish. 
However, tp+ belonging to D(Q + ) must also belong to 23(71"+) and since the functions in the 
form core for 7r+ behave as r a or nicer, we see that the term c^_ 1 r a ~ 1 e~ 1,e gets too singular 
to be in D(Q + ) if c+_ x ^ 0, and hence c^_ 1 must be zero. 

Similarly, a form core for t&l is given by e~ h (f>o, with <fio smooth. Functions in this form 
core have asymptotic behavior r~ a or r ~ a + 1 e l6 which forces c~ and c~_ 1 to be zero. 

Hence the parameters and uf 1 are given by Uq = c+/c^ a = oo, vf = cl_ a /c^_ l = oo, 
v o = c a/ c - a = and i/f = c{_ c Jc a _ 1 = oo. 

We conclude that the spin- up part of Ph can not be approximated in the sense of , while 
the spin-down part can. 
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